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Abstract 

The domain wall approach to lattice fermions employs an additional dimension, in 
which gauge fields are merely replicated, to separate the chiral components of a Dirac 
fermion. It is known that in the limit of infinite separation in this new dimension, 
domain wall fermions have exact zero modes, even for gauge fields which are not 
smooth. We explore the effects of finite extent in the fifth dimension on the zero 
modes for both smooth and non-smooth topological configurations and find that a 
fifth dimension of around ten sites is sufficient to clearly show zero mode effects. This 
small value for the extent of the fifth dimension indicates the practical utility of this 
technique for numerical simulations of QCD. 
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1 Introduction 



The anomalous breaking of the flavor singlet axial symmetry of QCD, ^7^(1), has important 
physical consequences. It is responsible for the relatively large mass of the r^' and 
effects the order of the finite temperature phase transition [Q. A central role in this is 
played by the index theorem which relates the winding of the gauge field with the 
number of zero modes of the Dirac operator. 

On the lattice the winding of the gauge field can not be defined unambiguously and the 
index of a finite-dimensional lattice Dirac operator is necessarily zero. Therefore, strictly 
speaking, the index theorem is not valid on the lattice. However, an approximate form 
of this theorem may still be present and govern anomalous effects. Some time ago this 
issue was investigated in for the case of staggered and Wilson fermions. A first step 
in that direction was to ask to what extent the lattice Dirac operator for staggered and 
Wilson fermions develops the appropriate number of zero modes for a fixed smooth gauge 
field background with given winding. It was found that for a particular background the 
appropriate number of zero modes was generated but they were not robust when high 
frequency noise was superimposed. Furthermore, the background chosen was of a very 
special nature and did not have the spatial variation present, for example, in an instanton- 
like background. Recent work 0, with more realistic instanton-like backgrounds, has 
demonstrated that staggered fermions do not develop zero modes unless the lattice spacing, 
a, becomes very small {a/D -C 0.1, where D is the instanton diameter). These difficulties 
suggest that for QCD, staggered fermions may fail to reproduce anomalous effects unless 
the lattice spacing is made quite small — considerably smaller than that used in present 
QCD thermodynamics studies, performed with Nt < 8 and staggered fermions. 

Current numerical results for the order of the finite temperature phase transition for two 
flavor QCD with light quarks show that it is not first order 0. This is consistent with the 
analysis in 0, provided Ua{^) is broken. However, current lattice discretizations manifestly 
break the symmetry of the classical action at finite lattice spacing and the corresponding 
discretized Dirac operators have difficulty seeing zero modes for non-smooth topological 
configurations. This raises the possibility that the present apparent second-order character 
of the two-flavor QCD phase transition may be a result of lattice artifacts rather than 
physical, anomalous symmetry breaking. A recent explicit calculation of an anomalous 
difference of susceptibilities near the two-flavor QCD phase transition, using a spectral 
sum-rule sensitive to zero modes, showed effects at or below the 15% level |p. However, 
in other work suggestions of anomalous zero-mode effects were seen. Clearly, a lattice 
fermion formulation whose action has the full global symmetry content of the continuum 
theory at finite lattice spacing would be a very useful tool for studying anomalous symmetry 
breaking effects. 



2 Domain Wall Fermions 

A novel approach with spectacular zero mode properties has been developed during the past 
few years. Domain Wall Fermions were first introduced in and a variant suitable for 
vector gauge theories was introduced in |]1T|, |T2|, |T^. (For reviews see and references 
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therein.) In this approach an extra space time dimension (henceforth to be called "s") is 
introduced with free boundary conditions on the two four-dimensional boundaries, s = 
and s = Ls — 1, where Lg is the extent of this new fifth dimension in lattice units. The 
fermion fields are defined in this extended space-time but the gauge fields are still defined 
on the ordinary space time and have no s dependence and no s component [jT5|. In this 
sense, the extra direction can be thought of as a sophisticated "flavor" space. Along the 
s direction the fermion field develops surface modes that are exponentially bound to the 
two free boundaries (domain walls) with the plus chirality component of the Dirac spinor 
on one boundary and the minus chirality on the other. As the size in lattice spacings 
(Ls) of this direction is increased the two chiral components get separated with only an 
exponentially small overlap remaining. For finite this overlap breaks chiral symmetry 
by an exponentially small amount and as Lg tends to infinity chiral symmetry is restored. 
Therefore, Lg provides a new parameter that can be used to control the regularization 
induced chiral symmetry breaking at any lattice spacing. For the first time the approach 
to the chiral limit has been separated from the approach to the continuum limit. 

An appealing aspect of the domain wall fermion formulation is the fact that the chirally 
symmetric, Lg —>■ oo limit can be analyzed in some detail using the overlap formalism 
P^ . Central to this overlap method are two Fock space states, designated \0h) and |0') 
in the notation of ||I2|, constructed from four-dimensional, single-particle, fermionic states. 
Because the gauge field is independent of s, it is possible to develop an s-independent 
transfer matrix (T) and associated Hamiltonian along this direction. Now in the Lg ^ oo 
limit, T^" becomes a projection operator onto the vacuum state, defined as \0h), the Fock- 
space state in which all the negative energy states of H are filled. The second state |0') is 
a much simpler, kinematic construction in which all single-particle eigenstates of position 
and 7^, with negative 7^ eigenvalue, are filled. 

A fermion Green's function for a given gauge field background can then be expressed as 
a simple matrix element of the appropriate number of creation and annihilation operators 
inserted between the states |0') and \0h)- In particular, the five-dimensional fermion deter- 
minant in the massless case is simply proportional to |(0'|Oii/) p which, for a finite space-time 
volume, can be calculated explicitly numerically. If the number of filled levels in \0h) and 
|0') is the same then |(0'|Oi^)p 7^ 0. However, if these filling levels differ, then |(0'|Oh)P = 
for zero mass and finite volume. This implies the presence of exact zero modes in the five- 
dimensional formulation. For a Green's function to be non-zero, an appropriate number of 
creation and annihilation operators must be inserted to balance the deficit. The deficit is 
naturally integer-valued and is defined as the index of the lattice Dirac operator. Therefore, 
this method provides a way to associate an index with the lattice fermion operator in the 
limit Lg ^ 00 but at fixed lattice spacing [O. 



More specifically it was found [15|, |16| that for classical backgrounds these zero modes 



are exact, that the deficit is equal to the winding of the gauge field and that they exhibit 
all the properties that are expected in the continuum. It was also found in that these 
modes are robust under the addition of high frequency noise. Using the overlap formalism, 
numerical simulations of the massless [|l^] and massive |]l8l vector Schwinger model gave the 
expected value for the anomalously generated t'Hooft vertex []T^. Furthermore, the index 
of the Dirac operator was calculated in SU(2) pure gauge theory slightly above the zero 
temperature crossover region |]T3. The index agreed within one sigma with the topological 
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charge as calculated in [20] indicating that the index theorem holds in a statistical sense. 
Similar studies were done in |2l| for pure SU(3) gauge theory. 

The above mentioned body of work indicates that the overlap (L^ oo limit of domain 
wall fermions) successfully incorporates exact zero modes at finite and relatively large 
lattice spacings. This suggests that with domain wall fermions, lattice QCD describes 
anomalous effects with regularization artifacts under firm control. Unfortunately, a direct 
implementation of the overlap formalism in QCD needs computing resources that are be- 
yond the capabilities of present day supercomputers. While a new proposal for reducing 



the computational cost of the overlap has been made |2^, [Q, it is not yet clear what the 



QCD computing requirement of this method will be. On the other hand, a straight-forward 
approach is to keep Lg finite. This method was used in to simulate the dynamical two 
flavor Schwinger model. In that work a detailed analysis of the Lg dependence indicated 
that for the massive theory the Lg = oo value of the chiral condensate and of the t'Hooft 
vertex was already reproduced within a few percent at Lg ^ 10. Furthermore, the rate of 
approach was consistent with exponential decay with a decay rate that became faster as 
the continuum limit was approached. Also, for an application to quenched QCD see [EH . 



3 Numerical Results 

These promising results indicate the important possibility of practical, chirally-consistent 
QCD simulations with domain wall fermions at finite Lg. In preparation for such simulations 
we wish to investigate to what extent the zero mode properties of the Lg = oo theory are 
maintained at numerically accessible values of Lg (~ 10 — 20). If a much larger Lg is needed 
then one would not be able to exploit these important features. On the other hand, if these 
properties are maintained at accessible values of Lg then anomalous effects can be studied 
with firm control over the finite Lg artifacts. 

Here we investigate this question using a classical instanton-like background [^ with 
prescribed winding of one unit. This is a compactified, singular-gauge instanton with origin 
in the center of a unit hypercube. The instanton field is cut off smoothly at some radius 
'"max so that it is entirely contained within the lattice volume. Specifically, we begin with: 

3 

A,{x) = -zJ:v""'^' , p(r)=Po(l-^)e(r^ax-r) (1) 

~{ + p{r)^ r^ax 

where is the gauge field potential, x^, is the space-time coordinate, r is the magnitude 
of X, X\ j = 1,2,3 are the first three Gell-Mann matrices, rj^'^'^ is as in |^, G is the 
usual Heavyside function and po is the instanton radius. Outside the fixed radius rmax the 
configuration is strictly a gauge transformation. This continuum field is then transcribed in 
the standard way to a lattice configuration of group elements, U^{x), defined on the lattice 
links X, IX. Finally, the lattice equivalent of the continuum transformation to singular gauge 
is applied: 

Uf,{x) ^ g{x)U^{x)g-\x + Uf,) , g{x) = ^^— (2) 

j=o m 
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where A" is the identity matrix and is a unit vector along the direction /i. Provided that 
the instanton center is not on a lattice site, this transformation is well defined everywhere. 
After the transformation, all links lying entirely outside rmax are equal to the unit matrix 
so the configuration exactly obeys the usual periodic boundary conditions. This instanton 
field is implemented on an 8^ lattice with po = 10, Tmax = 3 and on a 16^ lattice with 
pQ = 20, Tmax = 7. Thus, upou moviug from the smaller to the larger lattice, we have 
reduced the lattice spacing, measured in units of the instanton radius, by nearly a factor 
of two. In order to study the robustness of zero modes, random fiuctuations with a given 
amplitude ( are superimposed. In particular, at each link a different SU(3) matrix is 
constructed by exponentiating a linear combination of Gell-Mann matrices with random 
coefficients in the range [— C;C]- The gauge field at each link is then multiplied by these 
matrices. 

The operator used to study the effects of the winding of the gauge field background 
on the fermion sector is the chiral condensate (qq) calculated on that background. For 
staggered fermions (qq) = -^Tr[D^^ where D is the standard staggered Dirac operator 
and V is the four-dimensional volume. For domain wall fermions {qq) = j^Tr[{D~^)4^d] 
where D is the five- dimensional Dirac operator of [0], {D~^)4d is its inverse with the fifth 
dimension indices fixed so that it corresponds to a propagator between four-dimensional 



quark fields that are projections of the 5-dimensional fields as prescribed in [12|, and V is 



the four-dimensional volume. In particular, a four- dimensional Dirac spinor field is formed 
by combining the two upper spin components of the five-dimensional Dirac spinor field at 
s = with the lower two spin components at s = Ls — 1. This Dirac operator contains an 
explicit mass term that mixes the right and left chiralities with strength rrif. Antiperiodic 
boundary conditions along the time direction were implemented for both staggered and 
domain wall fermions. The inversion was done using the conjugate gradient algorithm. 
The stopping condition (ratio of the residual over the norm of the source) was set to: 10~^ 
for masses in the range [10~^,10~^], 10~^ for masses in the range (10~^,10~^], 10~^ for 
masses in the range (10~^, 10~% and 10~^ for masses in the range (10~^, 10~^]. 

The trace is over space-time, spin and color. A stochastic estimator was used to calculate 
the trace. To get reasonable estimates of the average and error one would have to use a 
large number (~ 50) of Gaussian random vectors. However, in this paper the interest is 
not so much in the actual value of the trace since it is only used as a device in studying 
the smallest eigenvalue (topological zero mode) of the domain wall fermion Dirac operator. 
To the extent that the Dirac propagator which enters (qq) is dominated by this single 
eigenvector of interest, the complete trace might be replaced by a single diagonal element 
taken in a random direction. This would give the desired trace, multiplied by the overlap 
between the random vector and the dominant zero mode. If the same random vector is used 
for all values of mo, rrif and Ls, it is expected that this overlap will be essentially constant 
and the variations seen will be those present in the full trace[6]. This strategy is followed 
in the calculation. However, in order to reduce the chance that the contribution of the 
single eigenvector of interest is accidentally suppressed by an unfortunate choice of random 
diagonal element, the same set of ten Gaussian random vectors is used to approximate the 
trace for all values of tjiq, mj and L^. The error bars shown in the figures represent the 
fiuctuations seen among these ten random vectors. Given this small number of samples, 
it is not certain that these error estimates are accurate. However, additional calculation 
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suggests that this particular sample probably underestimates the errors by less than a factor 
of two. 



At infinite Lg the explicit overlap formula for this formulation was derived in 
Using that it can be shown that {qq) diverges as (qq) ~ |z/|/m/, where u is the index of the 
Dirac operator as defined by the overlap formalism. At finite Lg the two chiralities do not 
completely decouple and therefore there is a residual mass. In free field theory this mass 
decreases exponentially with Lg and the effective quark mass is proportional to the sum of 
the residual mass and the explicit mass mj. In particular [1^ 



meS = mo{2-mQ)[mf + {l-mo)'] (3) 

where mo is the mass (domain wall height) of the five dimensional theory. In free field 
theory one flavor physics is obtained for mo in the interval [0,2) |TD[. In the interacting 
theory the boundaries of this interval will be renormalized. Assuming a similar modification 
of the effective quark mass for the case of the instanton-like background one would expect 
that at finite Lg, (qq) would behave as {qq) ~ |z^|/[mj + (1 — rriQ)^"], where uIq = mo(mo) 
is a "renormalized" domain wall height. 

As mentioned above, the index of the Dirac operator in the overlap formalism is nat- 
urally integer valued. A method to measure this index was developed in and used in 
p!9| , |T^, [^, [^. The index is half the difference of the number iV+ of positive minus 



the number of negative eigenvalues of the Hamiltonian H associated with the transfer 
matrix T. This number is the same as half the difference of the number of positive and 
negative eigenvalues of the operator D = •y^D^ where is the standard Wilson Dirac 



operator evaluated at a mass which is the negative of the domain wall height, i.e. —mo [|15 
||12|| . For mo < it can be shown that iV+ = N_ for any background gauge field. Therefore, 
by monitoring a few of the small positive and negative eigenvalues of D{mo) as mo is varied 
between zero and the positive value of interest, one can determine the number of positive 
eigenvalues that crossed zero and became negative and the number of negative eigenvalues 
that crossed zero and became positive. The difference of the number of the two types of 
crossings is the index. 

For the classical backgrounds that we studied we found that on the 8^ lattice the index 
changed from zero to one at mo ~ 0.28 and back from one to zero at mo ~ 2.14. For the 
16^ lattice the index changed at mo ~ 0.05 and mo ~ 2.01. ( Note this approach to the 
free-field values of and 2 is expected as one goes to the smaller lattice spacing implied 
by our 16^ instanton-like configuration. ) These values change by less than 1% when noise 
was added. Therefore, in these intervals one expects (qq) to diverge as rUes is made small. 
At the crossing points the transfer matrix has a unit eigenvalue and even at L<j = oo the 
two chiralities do not decouple |T^. Therefore, near a crossing larger values of Lg may be 
needed to see the expected 1/m behavior |jl2|. 

A "sketch" of the expected behavior of (qq) versus mj in the presence of a gauge field 
background with net winding is presented in Figure 1. In order to facilitate the comparison 
with staggered fermions, we will also use rrif for the usual staggered fermion mass. It is 
useful to analyze the mass dependence of (qq) by reference to the functional form of the 
spectral formula 

{qq) = [m/V]J2{l/i>^' + m')}. (4) 
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This continuum equation is exact for the case of staggered fermions and offers a useful 
framework for discussion of domain wall fermions. We can distinguish four distinct regions 
suggested by this functional form. In the large ruf region (I), we expect the propagator to 
be dominated by the mass term and therefore {qq) ~ ^/ruf. This behavior is expected and 
seen for both domain wall and staggered fermions. Next, we define region (III) for staggered 
fermions, as the mass range within which is small but larger than the smallest eigenvalue 
Ajnin. Here, Amin is presumably shifted away from zero by finite lattice spacing effects. For 
staggered fermions Amin is four-fold degenerate or near degenerate for the case where noise 
has been added. Region (III) is defined similarly for domain wall fermions, except the 
condition rrif > Amin is replaced by mj > mres, the residual mass due to the mixing of 
chiralities between the walls. (For the free field case, this is the (1 — mo)^'' term in equation 
^.) Although, domain wall fermions in the — >■ oo limit have Amin = its effects at finite 
Ls are cut off by this residual mass. Thus, for both staggered and domain wall fermions, 
one expects the small eigenvalue mode(s) to dominate the value of (qq) in region (III) and 
therefore {qq) ~ l/lrrifV) where V is the lattice volume. Region (II) is the crossover region 
between region (I) and region (III). In much of this region the mass is much larger than 
Amin but is Small enough to be relatively unimportant when compared with the rest of the 
eigenvalue spectrum. Therefore, referring to equation 4, one would expect {qq) ~ m/ in 
much of region (II). 

The expected behavior in region (IV) is different for staggered and domain wall fermions. 
For staggered fermions (solid line) Amin is small but not zero. Thus, in region (IV), where 
'^Z < Amin, the effects of the factor of in the numerator of eq. ^ will dominate and 
therefore one would expect {qq) mj. In contrast, for domain wall fermions (dashed 
line) and rnf < mres, we expect mj-es to play the role of m in eq. ^ so that in region (IV) 
{qq) ~ const. For domain wall fermions, region (IV) is dominated by the finite Lg, residual 
mixing of the two four-dimensional boundaries (see eq. |^). In the numerical results of 
this paper region (I) and part of region (II) are not present because the largest mass is 
rrif = 10~^. The focus is on region (III) where the divergent l/rrif behavior is expected 
and on the beginning of region (IV) whose onset signals the need for larger values of Lg. 

The numerical results are presented in Figures 2 — 8. In Figure 2 the staggered {qq) 
is plotted versus the quark mass m/ in the presence of the compactified, singular-gauge 
instanton background. Two lattice volumes and three different noise amplitudes, ( = 
0,0.01,0.1, are shown. For the case of the 8"^ lattice with no noise, one may be able to 
recognize divergent, 1/mj behavior in {qq) for 10~^ < rrif < 10~^ indicating the presence 
of a near-zero mode. This divergent behavior does not extend to smaller mj presumably 
because of the zero mode shift effect As can be seen when the lattice size is increased 
from 8^ to 16^ the instanton field becomes "smoother" (lattice spacing is reduced by a 
factor of two) and the divergent behavior becomes more pronounced now extending to the 
region 10~^ m-f < 10^^. For a detailed analysis the reader is referred to 0. However, 
when noise is added the divergent behavior begins to disappear. At noise amplitudes of 0.1 
the divergent behavior is not present in the 8^ lattice while in the 16^ lattice it has been 
significantly reduced. 

Figure 3 is the same as Figure 2 but now for the domain wall fermion {qq) with rriQ = 1.2. 
Results are shown for the two volumes 8'* and 16^ and for Lg = 4, 6, 8, 10. As can be seen, 
{qq) l/rrif for small mj. As expected, when rrif becomes smaller than some value, the 
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residual mixing between the two chiralities becomes the dominant contribution to meg and 
(gg) stops changing. This value is the border between regions (IV) and (III) sketched in 
Figure 1. For the 16^ lattice and for Lg = 6,8, 10 one finds that (qq) ~ l/mj for nif as 
small as ~ 10~^. For Lg = 4 there is no signal of a divergence because the effective quark 
mass is dominated by the finite Lg residual mass. For the 8^ lattice similar behavior is 
observed but now for Lg = 8, 10. For = 6 a divergence is observed but only for m/ down 
to ~ 10~^. When noise is added with amplitudes ( = 0.01 and C = 0.1 the behavior remains 
unaffected for all practical purposes indicating robustness under high frequency noise. For 
C = 1 the zero modes disappear (not shown here) but this level of noise is so large that it 
presumably destroys the winding. In particular, the index of the Dirac operator, as defined 
by the overlap formalism, was found to be zero. 

In Figure 4 the 16^ lattice of Figure 3 for = 0.1 and Lg = 10 is plotted again, this 
time with a (qq) = c_i/m/ fit for 10~^ ^ iTif < 10^^ superimposed. The coefficient c_i is 
1.09(7) X 10~^ and the per degree of freedom is ~ 0.5, demonstrating quantitatively the 
expected mj^ behavior. 

In order to study the mo dependence, (qq) versus mo is plotted in Figure 5 for fixed 
rrif = 5x 10~^. The vertical lines indicate the values of mo where 75!)^ (mo) has a crossing. 
One can see that near the crossing points larger values of Lg are required before (qq) 
becomes independent of Lg. However, there is a large range of mo in between the crossings 
for which (qq) does not change much when Lg is changed from 8 to 10. This means that no 
fine tuning of mo is needed even when large amounts of ultraviolet noise are added. The 
shape of the curves can be attributed to the wave function normalization factor which is a 
function of mo. In the free theory this factor is mo(2 — mo) ||18|| . 

To further verify robustness under changes of mo, (qq) versus m/ is studied for mo = 
0.75. As can be seen from Figure 5 this value of mo is in the onset of the region where 
the Lg dependence becomes stronger. The results for Lg = 4, 6, 8, 10, 12 are presented in 
Figure 6. Also, Lg = 24 is presented for the 8^ lattice. As can be seen, the zero mode 
effects are maintained and are robust under high frequency noise but larger Lg is needed 
before {qq) becomes independent of Lg. For the 16'^ lattice (qq) becomes independent of 
Lg for Lg = 12, but for the 8^ lattice this does not happen until Lg = 24. This is expected 
since the crossing point for the 8^ lattice is closer to mo = 0.75 than the crossing point for 
the 16^ lattice. 

If mo is allowed to get close to the crossing point one would expect that in order to 
maintain the zero mode effects much larger values of Lg may be needed. Figure 7 is the same 
as Figure 6, except that here mo is chosen a fixed distance of 0.2 from the crossing point. 
For the 8^ lattice mo = 0.5 and for the 16^ lattice mo = 0.25. The plots for Lg = 8, 12, 16 
and 24 are shown. Also for the 8"^ lattice the plot for Lg = 48 is presented. As can be seen 
for both volumes Lg = 24 is barely adequate to maintain the zero mode effect for masses 
5 X 10^^ < ^f- III order to maintain zero mode effects for masses 10"'* < my < 5 x 10"'^ 
an Lg = 48 is needed for the 8^ lattice. 

Finally, (qq) is plotted versus Lg in figure 8 for zero noise amplitude, fixed mj = 5 x 10~^ 
and for three values of mo = 0.5, 0.75 and 1.2 corresponding to circles, squares and crosses. 
The fits are to a function of the form l/[co + cie~'^^^'']. For the three values of mo the fitting 
range of Lg is [8,48], [6,48] and [4,48], the coefficient Cq = 0.183(8), 0.44(3) and 0.91(31), 
and the x^/dof is 4 x 10"^, 4 x 10"^ and 2 x 10""^. "(Note, these fits were performed by 
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minimizing the simplest which did not incorporate the strong correlation between the 
fluctuations at different values of L^. The presence of such correlations, resulting from using 
the same random source vectors for each L^, is reflected in the abnormally small values 
of x^.) Again, the quality of these fits demonstrates the expected l/rricQ ~ l/[mj + m^es\ 
dependence of (gg) but this time m/ is held fixed while m-^es ~ e'^^^-C-s jg varied. As expected, 
the decay rate is fast for trq = 1.2 but as mo gets closer to the value where has a 

crossing the decay rate becomes slower. Also, as expected from the discussion relating to 
equation 4 and figure 1 when mjcs becomes large the behavior changes from the monotonic 
behavior of regions IV and III to that of region II. For mo = 0.5 this can be seen for 
= 4, 6 and for mo = 0.75 it can be seen for Lg = 4. For mo = 1.2 presumably mres does 
not get large even for Lg = A and there is no change in the monotonic behavior. 



4 Conclusions 

In conclusion, for a classical, instanton-like background the domain wall fermion chiral 
condensate diverges as (qq) ~ 1/mj for rrif as small as 10^^ and Lg as small as ~ 10 lattice 
spacings. This behavior was observed on lattices of size 8^ and 16^ and was robust under 
random, high frequency noise with amplitudes as large as ( = 0.1. This l/rrif divergence 
was observed for a wide range of values of mo, indicating that there is no need for a 
fine tuning of mo. Furthermore, at fixed rrif and mo the chiral condensate approaches its 
Lg = 00 asymptotic value exponentially fast with a rate that becomes faster as mo is varied 
away from the point where the index of the Dirac operator changes. Therefore, domain 
wall fermions continue to show their spectacular Lg = 00 zero-mode properties even for 
values of m/ and Lg that are practically accessible to contemporary numerical simulations. 



Numerical simulations of QCD are currently under way p6[ to study the zero mode effects 



of domain wall fermions on gauge field backgrounds with realistic quantum fluctuations. 
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Figure 1. The expected functional form of (qq) vs. TJif. The sohd Une corresponds to 
staggered fermions and the dashed hne to domain wall fermions. 
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Figure 2. {qq) vs. for staggered fermions. The left columns are for an 8"^ volume and 
the right columns for a 16^ volume. The C = figures correspond to the compactified 
singular-gauge instanton background. The C = 0.01,0.1 figures correspond to the same 
background but random noise has been superimposed with amplitude 
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Figure 3. Same as Figure 2 but for domain wall fermions at mo = 1.2. Four different 
values of are shown with the circles, squares, crosses, and diamonds corresponding to 
= 4,6,8, 10. 
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Figure 6. Same as Figure 3 but for mo = 0.75. The circles, squares, crosses, diamonds 
and plus symbols correspond to Lg — 4, 6, 8, 10, 12 respectively. The star symbols on the 
8^ graph correspond to Lg — 24. 
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Figure 7. Same as Figure 6 but for uiq = 0.5 for tlie 8"^ and mo = 0.25 for tlie 16^. Tlie 
circles, squares, crosses and diamond symbols correspond to Lg — 8, 12, 16, 24 respectively. 
The plus symbols on the 8^ graph correspond to = 48. 
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Figure 8. (qq) vs. Lg for zero noise amplitude, fixed m/ = 5 x 10"'^ and for tliree values 
of mo = 0.5,0.75 and 1.2 corresponding to circles, squares and crosses. The fits are to a 
function of the form l/[co + Cie"'^^^*]. For the three values of itlq the fitting range of Lg is 
[8,48], [6,48] and [4,48], the coefficient Cq = 0.183(8), 0.44(3) and 0.91(31), and the xVdof 
is 4 X 10-2, 4 X 10-2 and 2 x 10-^ 
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